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We provide unambiguous theoretical evidence for the formation of correlation-induced isolated 
merons in rotationally-symmetric quantum dots. Our calculations rely on neither the lowest-Landau- 
level approximation, nor on the maximum-density-droplet approximation, nor on the existence of a 
spin-polarized state. For experimentally accessible system parameters, unbound merons condense in 
the ground state at magnetic fields as low as B* = 0.2 T and for as few as N = 3 confined fermions. 
The four-fold degenerate ground-state at B* corresponds to four orthogonal merons \QC) charac- 
terized by their topological chirality C and charge Q. This degeneracy is lifted by the Rashba and 
Dresselhaus spin-orbit interaction, which we include perturbatively, yielding spectroscopic accessi- 
bility to individual merons. We further derive a closed-form expression for the topological chirality 
in the form of a chiral spin current and use it to both characterize our states and predict the exis- 
tence of other topological textures in other regions of phase space, for example, at N = 5. Finally, 
we compare the spin textures of our numerically exact meron states to ansatz wave-functions of 
merons in quantum Hall droplets and find that the ansatz qualitatively describes the meron states. 
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I. INTRODUCTION 

Spin states of electrons confined in semiconductor 
quantum dots (QDs) are an appealing platform for quan- 
tum information science due to their long dephasing 
times and controllability 1 " 4 . States used for this imple- 
mentation must be long lived relative to the time-scales 
required to flip a qubit, and remain coherent despite sys- 
tem perturbations. Topological states are promising can- 
didates in that they are expected to have long decoher- 
ence times due to their global correlations, and be robust 
against local perturbation^. 

Skyrmions and merons are examples of topological spin 
textures that are pre dict ed to form in two-dimensional 
(2D) electron systems^'. Skyrmion excitations in bulk 
2D systems have been predicted in the v = 1 quan- 
tum Hall regime^^, and can condense into the ground 
state away from v = iP. Experimental evidence supports 
the existence of skyrmion spin textures in GaAs/AlGaAs 
quantum wella^Hll] 

A meron can be described as half a skyrmion; it con- 
tains a central spin oriented perpendicular to the 2D 
plane which transitions smoothly into an in-plane wind- 
ing away from the central spin. Unbound merons are 
not low-energy states in bulk 2D systems due to the 
prohibitive exchange energy of the in-plane winding. In 
finite-sized systems, such as the QDs we are considering 
in the present work, the textures are stabilized, and as 
we show below, they can condense in the ground state. 
Recent literature predicts the formation of merons in con- 
fined systems such as QDs in large magnetic fields^"^. 

In this work, we use configuration-interaction tech- 
niques to study a fully interacting QD system to pro- 
vide conclusive evidence for the formation of merons in 
the ground state. We find that merons form in systems 
containing as few as three particles^, at magnetic fields 



as low as 0.2 T, and away from the maximum density 
droplet regime. These states are degenerate. However, 
we show how this degeneracy can be lifted, and provide 
a method for predicting when merons will form in QD 
systems. 



II. CHIRALITY 

The chirality C n (C) of a vector field v(s), over a closed 
curve C in the direction of the unit vector n can be de- 
fined as 



C n (C) 



ds 



(v(s) x d s v(s)) 
\v±(s)\ 2 



(1) 



where v±(s) = n x v(s) is perpendicular to n and serves 
as a normalization factor; it is the global chiral character 
of the vector field in which we are chiefly interested. 

We confine ourselves to two-dimensional, rotationally- 
symmetric QDs centred at the origin, and we take the 
curve C to be a circle about the origin with radius r. 
Furthermore, we take our direction n to be perpendicular 
to the plane of the QD. Our vector field is the spin density 
S(r) = (S(r,9)). Taking the QD to lie in the x-y plane, 
Eq. ([I]) takes the form 
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(S{r,8) x d e S(r,0)) 
\S±(r,9)\ 2 



where S±(r,6) = zx S(r,0). 

We express our spin-density operator as 

S(r)=^t(r)<7 CTCT ^(r), 



(2) 



(3) 



where ffi a (r) is the field operator creating a fermion with 
spin a/2 (a = ±1) at position r and <x is the vector of 
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Pauli matrices. The chirality then takes the form of a 
chiral spin current 



with the chiral spin current density given by 



(4) 



-3 (S+dgS* + - S* + dgS + ) , (5) 



where S+ = (0?j (r) t/>^ (r)). As an example, consider the 

state |0) = (|£ 2 5 z )+ 7 |L; S^))/v^, where |L Z W S z (,) > is a 
general, correlated iV-particle state with orbital angular 
momentum (z- component) L z and spin angular momen- 



tum (z-component) S 
Eq. Q yields 



(') 



For these particular states 



C z (r) = (L' z -L z )(S z -S' z )6 ls ^\, 



(6) 



Perhaps not surprisingly, states of different orbital and 
spin angular momentum must be combined to produce a 
global chirality. C z is an integer winding number whose 
magnitude indicates the number of full 2-7T rotations along 
a closed curve C about the QD origin, and whose sign in- 
dicates the sense of rotation. In the interacting system, 
such windings can spontaneously occur at points of de- 
generacy in the spectrum. (Both spin and a component of 
orbital angular momentum are conserved in our system.) 
We show below that this can occur in parabolic QDs con- 
taining three strongly-interacting fermions at a point of 
four-fold degeneracy. We further show that the spin-orbit 
interaction splits this degeneracy into four spectroscopi- 
cally distinct merons \QC) with Q,C — ±1. 



III. THEORY 

Our system consists of N interacting fermions of charge 
e, bound to a 2D plane and laterally confined by a 
parabolic potential. The 2D Hamiltonian used to de- 
scribe this "standard model" is 



N 



N 



1 

2^e|fi 



H 



so, 



(7a) 



where e is the dielectric constant of the medium and h is 
the single-particle Hamiltonian describing harmonic con- 
finement in a perpendicular magnetic field: 



1 1 

- — -P (r) + -m wzr 
2m 2 



2 -2 



qHbS z B z , 



(7b) 



where r = (£, y) is the position operator, m* the effec- 
tive mass, and ujq the parabolic confinement frequency. 
We work in a finite magnetic field perpendicular to 
the plane of the QD, and so P = p + -A(r), with 
A (f) = |r (— y,x, 0) in the symmetric gauge. The fi- 
nal term in Eq. ( 7b ) is the Zeeman energy and plays no 
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FIG. 1: Low-lying spectra for N = 3 interacting particles 
with 2D harmonic confinement, Eq. |7f. L z and S quantum 
numbers are displayed for each state. The field B* at which 
merons form in the ground state is marked. 



significant role in what follows due to the relatively small 
magnitude of both the <?-factor (g sw —0.44 in GaAs) and 



the magnetic field (B < 1 T). The final term in Eq. (7a 



is the spin orbit interaction, which we deal with pertur- 
batively in a subsequent section. In what follows, we fix 
the material parameters in Eq. ([T]) to typical values for 
GaAs heterostructures: m* = 0.067m e , e = 12.4, and 
hiu a = 1 me\fiM. Finally, we take ft/2 = 1. 

Neglecting for now the spin-orbit interaction, we ob- 
tain the correlated eigenstates and the energies ofEq. Q 
by configuration-interaction methods. The single particle 
system, Eq. ( 7b ) , can be solved exactly yielding the well 
known Fock-Darwin spectrum^ with states \mns) la- 
belled by harmonic oscillator numbers m, n = 0, 1, 2, . . ., 
and spin s = ±1/2. We use these states to build up 
a basis of antisymmetrized many-body states, and write 
the Hamiltonian as a block-diagonal matrix with blocks 
segregated by the conserved quantities of particle num- 
ber N, total spin S, z-component of spin S z , and orbital 
angular momentum L z . 

The simplest system exhibiting meron textures occurs 
already for N = 3 confined particles. The low-lying spec- 
trum is shown in Fig. [T] For the system parameters given 
above, a ground state degeneracy exists at a field of B* = 
0.365 T. This degenerate manifold is spanned by the four 
spin- 1/2 states \L Z ,S Z ) = | — 1, ±1/2), |— 2, ±1/2). In the 
harmonic oscillator basis, these correlated states typically 
contain several thousand Hartree-Fock states in order to 
reach convergence in the energies to within 0.05%. 



IV. MERON TEXTURES 

The degenerate ground-state subspace at B* described 
above contains merons. An explicit form which is impor- 



FIG. 2: (Color online) Two of the four meron textures for 
the three-particle states given in Eq. J8j. The spin densities 
(Sx(r)) and (S y (r)) are represented by the vector field, while 
(S z (r)) is represented by the color bar. The remaining two 
configurations are obtained from these by a local in-plane spin 
rotation of it. The textures shown were obtained at the four- 
fold degeneracy point B = B* (see Fig. [TJ. 



tant in what follows is 
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FIG. 3: (Color online) Comparison of spin density compo- 
nents for meron \QC) = I 1 1) of Eq. (|8| (solid lines) and the 
semiclassical result of I!cf. " (dashed lines) for 6 — 0. The 
S x (r) curves vanish at the origin, and S y (r) is zero in both 
cases here. Inset: Particle density for the |1 1) meron (solid) 
and semiclassical (dashed) result. 
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These four states for C = ±1,Q = ±1 are orthogonal and 
yield the spin textures (S{r)) shown in Fig.jij These spin 
textures are plotted using the real-space wave functions 
of the 2D harmonic t raj {20. Previous worlfl^E^I has fo- 
cused on lowest-Landau-level physics or on semiclassical 
approximations. The present numerically exact results 
show that merons can exist far beyond the semiclassical 
regime; right down to the extreme quantum limit of very 
few confined particles where correlations are strongest. 
In particular, the three-particle ground states shown in 
Fig. [2] are not spin-polarized states, nor do they corre- 
spond to the v — 1 maximum density droplet. The states 
in Fig. [2] correspond to a filling factor of 2 > v > 1. (The 
angular momentum of the ./V-particle maximum density 
droplet is L z = —N(N - l)/2.) 

A semiclassical (SC) ansatz for a single meron in a 
quantum Hall droplet is given in Refjl™ in terms of 
spinors and Laughlin wave functions. Figure [3] shows 
a direct comparison between the spin textures {Sir)) 
for the \QC) = |11) meron presented here and the 
SC ansatz for N = 3. Results are shown in units of 
£ = WK/Zrrt*LJ, where u = ^uJq + and u c is the 

cyclotron frequency. The two chief distinctions between 
the curves are the relative suppression of both charge 
(primarily at the origin) and spin polarization through- 
out the dot in the present work. Both effects can be 
attributed to the strong correlation effects captured in 
the present treatment. The suppression of charge at the 
origin is a result of the greater repulsion among the par- 
ticles. The correlation-induced mixing among numerous 
Landau-levels and angular-momentum and spin-resolved 
Slater determinants likewise has the effect of reducing the 
spin polarization throughout the dot. Indeed we find a 
conspicuous lack of full polarization throughout the dot. 



V. SPIN-ORBIT COUPLING 

To linear order, the spin-orbit interaction splits the de- 
generacy among the four merons states \QC). Remark- 
ably, we find the Rashba term couples only to merons 
with positive chirality, whereas the Dresselhaus term cou- 
ples only to those with negative chirality. Explicitly, the 
linearized spin-orbit interaction may be written as EMU] 

Hso = (3 {-o- x P x + a y P y ) + a (a x P y - cr y P x ) , (9) 

where a and (3 are the Rashba and Dresselhaus spin- 
orbit coupling strengths, respectively. Within the lowest- 
energy degenerate subspace, the spin-orbit-induced en- 
ergy splittings Eqc associated with each state \QC) are, 

E QC = S ^[(1 + C)a+(1-C)p], (10) 

where, 

£c = A<-2, f I [(1 - l)S c P (1 + l)S c a] I - 1, -f ), 

with A = l/(V2^o) an d Q,C = ±. Here, a and w in- 
dependently lower the orbital angular momentum of the 
\L Z , S z ) states, while S± raises or lowers the spin. Also, 
v = \/l + 4:UJq /lj^ . Note if either a or (3 is zero, or if 
a. = {3, a two-fold degeneracy remains. Otherwise, the 
degeneracy is completely lifted. Since the Rashba term 
is to some extent tunable through externally-applied gate 
voltages^2Hl2l j Eq ( 10 1 demonstrates a measure of exper- 



imental control over the merons. 



VI. PREDICTIONS 

Equation ^ predicts that degenerate manifolds with 
states of different orbital and spin angular momentum 
contain quasi-topological winding spin textures. Further- 
more, the magnitude of the winding number associated 
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FIG. 4: (Color online) Low-lying spectra for TV = 5 inter- 
acting particles with 2D harmonic confinement, Eq. |7J. The 
quantum numbers L z and S of selected levels are shown. B* 
marks the field at which merons form in the ground state. 



where L z £ L' z , \S Z - S' z \ = I, c = C/\C\, and q = Q/\Q\. 
The existence criteria for the quasi-topological wind- 
ing textures described by Eq. ( |l~2"| ) is satisfied through- 
out the phase space of the QD. For example, Fig. [4] 
displays the lowest-lying spectrum for the five-particle 
system using the same experimental parameters listed 
above. The four-fold degeneracy at B* = 0.16 T oc- 
curs between the states \L Z ,S,S Z ) = | — 1,1/2, ±1/2) 
and | - 4, 1/2, ±1/2). According to Eq. the quasi- 
topological winding states which exist at this point have 
a winding of three. Indeed this can be seen in Fig. [5] 
where the individual spin components for a meron state 
with QC quantum numbers \QC) = |1 — 3) are plotted. 
The large, off-centre peak in the (S z (r)) distribution is 
due to the increased Coulomb interaction strength rela- 
tive to the three-particle system. 



with these states is equal to the difference between the 
orbital angular momentum of the degenerate states. In 
the case of a four-fold degeneracy, Eq. ^ generalizes to 

\QC) = ^(\L z ,S z )-iq\L' z ,S' z )) 

+ ]±5(\L z ,S z } + q\L z ,S z }), (12) 
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